Cavity QED treatment of scattering-induced efficient free-space excitation and 
collection in high-Q whispering-gallery microcavities 
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Whispering-gallery micro-cavity laser possesses ultralow threshold, whereas convenient free-space 
optical excitation and collection suffer from low efficiencies due to its rotational symmetry. Here 
we analytically study a three-dimensional microsphere coupled to a nano-sized scatterer in the 
framework of quantum optics. It is found that the scatterer is capable of coupling light in and out 
of the whispering-gallery modes (WGMs) without seriously degrading their high-Q properties, while 
the microsphere itself plays the role of a lens to focus the input beam on the scatterer and vice versa. 
Our analytical results show that (1) the high-Q WGMs can be excited in free space, and (2) over 
50% of the microcavity laser emission can be collected within less than 1°. This coupling system 
holds great potential for low threshold microlasers free of external couplers. 

PACS numbers: 42. 60. Da, 42.55.Sa, 42.50.Ct 



I. INTRODUCTION 



Whispering-gallery mode (WGM) microcavities repre- 
sent one of the most promising candidates for a wide 
range of fundamental studies and applications, includ- 
ing cavity quantum electrodynamics, cavity optomechan- 
ics, microlasers, filters and biological sensors (for reviews, 
see Unfortunately, due to the rotational symme- 

try, they suffer from inefficient coupling with the out- 
side modes, which limits their applications, especially 
in microlasers One of the solutions is the tapered 
fiber coupling method (fHH , which possesses nearly unity 
efficiency. Nevertheless, convenient free-space coupling 
without near-field couplers is e agerly required because 
of the experimental limitations [9Hll| . For example, the 
external couplers are not convenient at low temperature 
chambers; for a higher-index-material resonator (l2l. fl3j. 
its coupling with the tapered fiber is inefficient due to 
the phase mismatch. Alternatively, deformed cavities 
(also named as asymmetric resonant cavity, ARC) are 
proposed fl4rl30| , because they allow high- efficiency free- 
space excitation and directional emission. Latest devel- 
opments in ARC studies include Limacon-shaped cavity 
[3ll434| and circular disk cavity with a linear |35[ or point 
defect in it [H,[33] , but the emission divergence angles are 
still too broad. Very recently, highly directional outputs 
are obtained in an elliptical microdisk with a notch at 
the boundary (3|| , and by placing a nanoparticle into the 
evanescent wave region of microcavities (39|. These in- 
vestigations focused on two-dimensional (2D) microcav- 
ity systems by resorting to numerical simulations. 
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However, 2D microcavities exhibit relatively low qual- 
ity (Q) factors in experiment and there are significant en- 
ergy losses in the perpendicular dimension. Thus, there is 
an interest to employ 3D microcavities, whereas it is dif- 
ficult to perform numerical simulations. Here we present 
a cavity quantum electrodynamics (QED) treatment of 
3D microsphcre-scatterer coupling system, and derive an- 
alytical expression of the free-space excitation and the 
emission directionality. Our results explicitly reveal the 
underlying physics and is also suitable for 2D case. 

The paper is organized as follows. In Sec II, we briefly 
describe the microsphere-scatterer coupling system. In 
Sec. Ill and IV, we investigate the free-space excitation 
of WGMs and free-space collection of scattered lasing 
modes, respectively. Conclusions are presented in Sec. 
IV. 



II. SCATTERER-MICROSPHERE SYSTEM 

Figure [IJa) illustrates a schematic of the present sys- 
tem. A spherical subwavclcngth scatterer (radius r s ) lo- 
cates on the surface of a microsphere (radius R) which is 
doped with gain medium for the microlascr applications. 
A Gaussian pump laser beam (vacuum wavelength Ai, 
satisfying r s < Ai < R) with the polarization in x axis 
direction propagates along —z axis, and is incident to 
the microsphere. Here we have established the Cartesian 
and spherical coordinate systems with the scatterer lo- 
cated at the origin O, as sketched in the bottom left of 
Fig. [Ua). The subwavelength scatterer can be treated 
as a dipole [4(| EH , with the dipole moment induced by 
the electric fields of the input Gaussian modes, the exci- 
tation and lasing WGMs and the reservoir modes in the 
free space. The Rayleigh scattering results in the inter- 
action among these modes, by which the input photons 
are scattered into the excitation WGMs, and the lasing 
photons in WGMs are scattered into the reservoir modes. 
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FIG. 1: (Color online) (a) Schematic illustration of the 
scatterer-microsphere coupling system (not to scale). Bottom 
left: the corresponding Cartesian and spherical coordinate 
systems, with 9 the emitting angle and <f> the azimuth an- 
gle, (b) Finite element simulations (COMSOL Multiphysics) 
of the directional emission mode. The color represents the 
square root of the electric field's absolute value l^l 1 / 2 (For 
clarity, we use |_E| 1/ ' 2 in order to increase the contrast). The 
red arrows indicate the direction of power flow. In the simula- 
tion we use the radii R = 10 /im, r s = 200 nm, microsphere's 
refractive index n = 1.7. 



The scattered lasing photons are collimated by the mi- 
crosphere, giving rise to directional emission. In Fig. 
QJb) the simulation results of wave approach is presented, 
demonstrating this directional emission mode. As a uni- 
versal paradigmatic approach, the quantum treatment 
of Rayleigh scattering is widely used and well demon- 
strated [io, I42U45I . This quantum treatment is necessary 
because in such a scatterer-microsphere coupling system 
there exists a phenomenon similar with the Purcell ef- 
fect, since the high state density of WGMs causes the 
enhanced scattering between the WGMs and the free- 
space reservoir modes [4(| EH . In the following we adopt 
this quantum approach to analyze the free-space excita- 
tion of WGMs and the free-space collection of scattered 
lasing modes, respectively. 



III. FREE-SPACE EXCITATION 

In this section, we develop the theoretical model to 
describe the scattering-induced coupling among the sur- 
rounding optical modes. By writing down the total 
Hamiltonian, we derive the equations of motion (quan- 
tum Langevin equations) for the free-space excitation 
process. Then we find analytical expressions for the cou- 
pling coefficients and the excitation efficiency. We also 



address the focusing effect of the microsphere, which en- 
sures efficient free-space excitation. At last specific ex- 
amples are presented. 



A. General Hamiltonian and Quantum Langevin 
equations 

Under the rotating wave approximation, the total 
Hamiltonian for the free-space excitation process can be 
written as 



rrtot rxf , rrf , rrf 

11 1 — fl cav ' a in ' -"res 

' cav-cav ' cav— in ' c 



cav— res ' 



(1) 



Here the superscript f (i) labels free (interaction) terms. 
The first three terms (h = 1) 
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(3) 
(4) 



describes the free radiation parts, where a m , b(to) and 
Cj denote annihilation operators of the m-th excitation 
WGM, the input modes and the j-th reservoir mode, 
with frequencies wi. m , to and loj, the commutation re- 
lation [a m ,a} m ,] = 5 m , m >, [b(u),tf(u')] = 5(u> - w'), 

[cj,ct] = Sjji, respectively. For the WGMs in Eq. (0), 
we can neglect high order WGMs and focus on the fun- 
damental WGMs, since they distribute around the equa- 
tor of the microsphere, possess the minimum mode vol- 
ume, the maximum Q factor, and are the typical lasing 
modes in actual experiments. Also, In the system de- 
scribed in Fig. HJa), only transverse electric (TE) WGMs 
can be efficiently excited, and transverse magnetic (TM) 
WGMs can be safely neglected since the electric field of 
the TM WGMs are almost orthogonal to that of the in- 
put modes. Therefore, in Eq. @ the summation index 
to runs through clockwise (CW) and counterclockwise 
(CCW) propagating fundamental WGMs, with a degen- 
erate frequency denoted as cji jTO = U)\. 

The interaction terms of Eq. ([1} are given by 



H 1 

cav— cav 
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(5) 



m 



cav — m 



/ [igi„, m (w)6 t (w)a m + H.c] du, (6) 
^2(9i, m .,Mn c j +H.c). (7) 



TJl 
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Here if' describes the interaction between the cav- 

ity modes, which results in the scattering into the same 
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(to = to') or the covmtcrpropagating (m ^ to') WGMs 
with amplitude coupling strengths gi. m ,m r , H 1 cav _ in 
represents the scattering between the input modes 
and the WGMs with coefficients g injTO (w); H 1 cav _ ICS 
describes WGMs-reservoir scattering with coefficients 
9i,m,j- Without loss of generality, in the following the 
coupling coefficients <?i, m ,m', 5in,m(^) and gi, m ,j are as- 
sumed to be real for notational convenience. The de- 
generate CW and CCW WGMs have the same coupling 
strength in the scattering process, and thereby we denote 

9i, m ,m' = 9m,m(u) = 5inM and gi, m ,j = 9i,j- Note 
that in Eq. ([1]) we have neglected interacting terms like 

H L-in> ^res-res) H L-res' sinCe the y do llot directly affect 

the cavity modes, and have minor effect on the system 
dynamics. 

The cigenmodes of the system are the symmetric and 
antisymmetric standing modes, given by a± = (acw ± 
accw)/v / 2. Using the Heisenberg equations and the 
Markov approximation [46l l47j , we obtain (see Appendix 
A for details) 



^± = -i(uJi - 2g 1 )a + - (n in + K R )a + 



where 



— y/2n in b hl - £, 

K in = 27rft 2 n (wi - 2g 1 ), 

kr = 27r^ - oji + 2gi) 



(8) 

(9) 
(10) 



denotes the input- WGMs energy coupling strength and 
the decay of the WGMs induced by the Raylcigh scat- 
tering to the reservoir. b- ln is the input field and £ is the 
noise operator relate to the reservoir. 



B. Scattering coefficients and excitation efficiency 



for a microsphere; /i(R) = |-E , i(R)/£ , i jmax | is the nor- 
malized field distribution function of the WGMs. The 
quantized electric field of the input beam at position 
(r,z) is given by E in (r,z) = [^ n +) (r, z) + E£\t, z)]e x , 
where r = (xe x , ye y ). The positive frequency component 
reads 



2e U ; 



7in,fe(r)6 fe e 4 



(14) 



where Sb denotes the relative permittivity of the sur- 
rounding medium, bk is the annihilation operator of the 
&;-th mode, Vk (/in,fc(r)) is the corresponding mode vol- 
ume (field distribution function) . The expression can be 
rewritten as an integral form [47| 



f in (v)b(uj)e^ z . (15) 



Here c is the speed of light in vacuum; A(z) is the mode 
area on x — y plane, given by 



: / g (r)|£ in (r)| 2 rfr 2 
max[e(r) \E in (r)\ 



A(z) 



(16) 



For Gaussian beam A(z) = JJ^ e" 2 ^ +v dxdy = 

irw(z) 2 /2, where w(z) is the spot radius at z. The quan- 
tized electric field of the reservoir is given by Ej(R) = 

[£j +) (R) + £H(R)]e i5 where 



K (+) (R) 



' Huij 



(17) 



In the full quantum theory, gi n (ui) and g\ can be cal- 
culated as follows. The quantized electric field of the 
excitation WGMs at position R is given by Ex(R) = 

[e[ + \h) + e[ \ll)]e x . Here e x is the unit vector along 
the x-axis direction, 



E[ + \R) 



/i(R)(acwe 



tki-R 



«ccwe ) 



(11) 

is the positive frequency component of the field and 
^"'(R) is its adjoint; ki is the wave vector of the CW 
mode; Eq is the dielectric permittivity of the vacuum and 
e c denotes the relative permittivity of the microsphere; 

Je(R)|£i(R)| 2 rfR 3 



max[e(R) |£i(R)| ] 



(12) 



is the mode volume of the WGMs, which can be calcu- 
lated as m 



\2Tm) 



(13) 



with Vj , kj and being the mode volume, wave vector 
and the unit vector along polarization direction of the 
j-th reservoir mode. 

The interaction via scattering yields the Hamiltonian 

MM M M 



where 



H i = ~p B -B a , 



E B =Ei(0)+E in (0,0) + Ej(0), 
p s = e aE s 



(18) 



(19) 
(20) 



are the total electric field at the position of the scatterer 
(the origin O of the coordinate system) and the total 
polarization of the scatterer; a = 4-7rr 3 (e s — l)/(e s + 2) is 
the polarizability of the spherical scatterer with £ s being 
its relative permittivity. Note that for Gaussian beam 
the maximum electric field is at the center of the beam, 
yielding /i n (0) = 1. Using the above equations, we obtain 
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the coupling coefficients as 



9i 



91,3 



i a/l(0) V^Sf^ 



(e 



(21) 
(22) 
(23) 



and thereby the in-coupling strength can be obtained as 
a 2 /i 2 (0Vi(wi-2 5l ) 



KB 



4e c cVi4. 
(n 5 + l)aV(^ 1 -2g 1 ) 3 / 1 2 (0) 
127rc 3 e c Vi 



(24) 
(25) 



where A s — with w s being the spot radius at 

z = plane (where the scatterer is located at). 

The excitation efficiency can be defined as rj = 1— T m i n , 
where T m ; n is the minimum value of the transmission. 
From Eq. <j8j> , we obtain (see Appendix B for details) 



n 



4Atin(^0 + «r) 
(«0 + Kin + Kr) 2 



(26) 



C. Focusing effect of the microsphere and specific 
examples 

To realize efficient coupling, the input beam should 
be focused into a small spot on the scatterer, since K- ln 
is in inverse proportion to the mode area A s , as shown 
in Eq. (|2~4")) . In fact, the microsphere is a natural op- 
tical lens which possesses ultrashort focal length, capa- 
ble of focusing the light spot significantly. Numerical 
simulation shows that circular dielectric cylinders illumi- 
nated by a plane wave can generate nanojets with waists 
smaller than the diffraction limit [4!| [50| ■ Here we use 
Gaussian beam input and analytically treat the problem 
using Gaussian beam transform laws. The sphere can be 
viewed as a thick lens, with the focal length 



nR 



2(n-l)' 



(27) 



where n is the relative refractive index between the sphere 
and the environment (n = y/E^). As depicted in Fig. 
Ufa) , the input light is assumed to be Gaussian beam with 
a waist radius wo (wq < i?), and the distance between the 
beam waist and the center of the microsphere is s. After 
being focused by the sphere, the beam waist becomes 



Fw 



(28) 



with its distance to the center of the microsphere given 
by 



(s-F) 



f)+z If 



(29) 




(b) 



Input beam waist w (\xm) 



FIG. 2: (Color online) (a) Illustration of the microsphere- 
induced focusing effect, with wo being the input beam waist 
and w s being the resulting spot radius at z = plane, after 
being focused by the microsphere, (b) w a vs. wo for differ- 
ent cavity refractive index n. Inset: spot radius minification 
w s /wo as a function of wq. The dotted vertical lines indicate 
too = 5 /im. Here the microsphere's radius R — 10 /J,m. 



where zr, = ttivqu/X is the Raylcigh range. Then the 



spot radius at z = plane is 




(30) 



with z^ = TTw'^n/X. 

In Fig. Hfb) we give a specific example, where R = 10 
/im, s = 0; for dopant Er 3+ , the excitation wavelength 
Ai = 977 nm (rasing wavelength A2 = 1550 nm). We 
plot the resulting spot radius w s as a function of the 
input beam waist wo for different refractive index n = 
1.5, 1.7, 1.9. It is shown that for relatively large input 
beam waist, the microsphere is able to focus the beam 
intensely, especially when n approaches 2. For F zr, i. 
e., Wq 3> Ai?/[27r(n — 1)], we obtain w s ~ (|2 — n\ /n)wo, 
which indicates a \2 — n\ /n times decrease of the sport 
radius, as further plotted in the inset of Fig. [2Jb). This 
is consistent with the ray optics predictions. Therefore, 
it is of great advantage to make use of the cavity itself as 
a micro-sized lens. 

Now we study the in-coupling strength re; n and 
Rayleigh scattering induced decay kr as a function of the 
scatterer's radius r s as shown in Fig. [Ha)-(b), where we 
have set wq = 5 ^m, and other parameters: e s = 12 (sil- 
icon), /i(0) = 0.4 (We will use there parameters unless 
specified). The free-space excitation efficiency for differ- 
ent radii of both the scatterer and the microsphere are 
presented in Fig. [2Ic)-(d), which shows more than 10% 
excitation efficiency can be obtained for suitable parame- 
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FIG. 3: (Color online) (a)-(b) In-coupling strength K[ n and 
Rayleigh scattering induced decay kr as a function of the 
scatterer's radius r a for different cavity refractive index n. (c)- 
(d) Excitation efficiency rj as a function of r s and _R. In (a)- 
(c), we use R = 10 ^m, corresponding to the dotted vertical 
line in (d); in (d), r s — 50 /im, corresponding to the dotted 
vertical line in (a)-(c). 



ters. In Fig. Etc), for small scatterer, the cavity intrinsic 
decay kq dominates over n ln and kr, which results in low 
excitation efficiency. For large scatterer, kq can be ne- 
glected compared with kr, yielding a constant excitation 
efficiency decided by K; n and kr. In Fig. |3|d), smaller 
microspheres possess smaller mode volumes, resulting in 
larger K- m and kr, and thereby larger excitation coeffi- 
cient. In addition, smaller microspheres have stronger 
focusing effects, leading to larger « in . To obtain a better 
excitation efficiency, we should increase K; n and mean- 
while decrease kq and kr. This can be realized by using 
smaller input beam waist and by using microcavity with 
proper refractive index, which lead to small mode area 
A s , as discussed above. 



IV. FREE-SPACE COLLECTION 
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FIG. 4: (Color online) Out-coupling strength n out as a func- 
tions of the scatterer's radius r a (a) and microsphere's radius 
R (b) for different cavity refractive index n. 



coefficients in the above section, we obtain 



(n 5 + l)a'u; 2 (u J2 -2g 2 yfi(0) 
12tt c 3 e c V 2 



(31) 
(32) 



where 



92 



^2/|(0) 

2e c V 2 ' 



92, J 



_i a/a (o)./J^z_(g !B .e,-) J 



(33) 
(34) 



For free-space collection process, the Hamiltonian is 
similar with that of free-space excitation process (Eq. 
([1])) by dropping the terms containing the input modes. 
Mention that the energy scattered from the lasing WGMs 
into the reservoir modes is just the output of laser emis- 
sion. Quite different from the ordinary case in which this 
kind of scattering leads to pure damping and is always 
harmful, here it is a kind of useful resource and plays 
a key role in obtaining directional laser emission. This 
scattering offers an interface between the WGMs inside 
the microcavity and the optical modes outside the cav- 
ity, where we use the coefficient K out to denote the out- 
coupling strength. Following the calculation of scattering 



and V 2 (/2(f))) is the mode volume (field distribution 
function) of the lasing WGMs. 

Figure 0] plots K out as a function of the radii of the 
scatterer and the microsphere. For r s = 50 nm and R = 
10 /im, K out is several mega Hertz, which corresponds to 
Qout = W2/« out ~ Therefore, the high-Q properties 
of the lasing modes can be maintained. 

In the following we analyze the emission directional- 
ity originating from Rayleigh scattering and collimating 
effect of the microsphere. Finally the emission direction- 
ality and energy collection ratio for various parameters 
are present. 




FIG. 5: (Color online) (a) Spherical plot of u{6, <j>). (b) u(9, 4>) 
vs. 8 for various (f>. Inset: The corresponding polar plot of 
u(0, 0). Here we use n = 1.5. 



A. Scattering directionality 



FIG. 6: (Color online) (a) Illustration of the microsphere- 
induced collimating effect, (b) vs. 6 for n = 1.5, 1.7, 
1.9. The dotted vertical lines correspond to the critical angles 
of total reflection for n = 1.5, 1.7, 1.9 (from right to left). 
Inset: Transmission as a function of 9 for <j> — 0° (solid), 90° 
(dashed), here n = 1.5. 



The factor • e\,- in Eq. (|34[) indicates that 
<?2j depends on the direction, resulting the direction- 
dependent out-coupling coefficient K,(9,<f>), which satis- 
fies K out =// K(0,4>)dn, with fl being the solid angle 
and dCl = s'm9d9d(f>. After normalization, we can define 
u(9, (j)) = k(6, <f>)/K ou t, which can be calculated as 



«(M) 



3ra a (l-sin^ 6>(W 

47r(n 5 + l) 
3(1— sin" cos (j)) 
. 47r(n 5 + l) 



, 0° < 6 < 180°, 
180° < 9 < 360°. 



(35) 



The normalized out-coupling coefficient u{6, 4>) repre- 
sents the angular distribution of output energy. In Fig. 
[5fa) we present u(9,4>) for given azimuth angle = 0°, 
45°, 90°, respectively. Note that for 0° < 6 < 180° the 
environment is dielectric cavity with permittivity of e c , 
while for 180° < 9 < 360° that is vacuum with permittiv- 
ity of 1. As shown in Fig. [5j the light trends to be scat- 
tered to z axis (9 = 0°,180°), but the scattering along 
+z axis is much stronger due to the asymmetry environ- 
ment. Note that for 4> = 90°, the scattering is uniform 
in the same environment, since in this case the scattered 
light has the same polarization with the WGMs. 



B. Collimating effect of the microsphere 

Although the scattered light trends to propagate along 
z axis, the divergency angle is too wide. To obtain bet- 
ter directionality, once again we can make use of the mi- 
crosphere itself, which behaves as a thick lens. As de- 
picted in the inset of Fig. HUb), the emitted light for 
0°< 9 < (180°/7r)arcsin(l/n) (the critical angles of to- 
tal reflection) passes trough the microsphere, and finally 



yields the output angle 0, given by 

© = f{9) = \26 - arcsm(nsm0)| 



(36) 



Note that for (180%r) arcsin(l/n) < 9 < 90°, the light is 
totally reflected on the microsphere surface. From Fig. 
[6jb) we can see that the output angle O has a much small 
divergence than the emitting angle 9, which stems from 
the collimating effect of the microsphere. 

The output energy density function can be obtained as 



p(6,0) = ^r(r 1 (e),0Hr 1 (e), 

where 



1 de 



T(e,<t>) = T p (6,< 



(37) 



(38) 



T p (T s ) is the transmission for p-polarization (s- 
polarization) component, calculated from Fresnel for- 
mula; 9 = / _1 (0) is the inverse function of f(9), which 
is a multiple valued function, and means the summa- 
tion over each section of the multiple valued function. 
Note that p(Q, <j>) has singularities, thus the full-width- 
of-half-maximum definition of the divergence angle fails. 
To quantify the emission directionality, we define a half- 
energy angle Q 1/2 , given by P(Q 1 / 2 ,4>) = 1/2, where 



P(0,0)= / p(0'»de' 

'0 



(39) 



is the energy ratio (energy distribution function), repre- 
senting how much energy distributes in the interval [0,0]. 
This half-energy angle represents that the output angle 
of half light is smaller than 0i/2- 
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As two specific cases, the energy density function 
p(0,</>) and the energy ratio P(0,</>) for n = 1.5 and 
1.9 are plotted in Fig. Elb)-(c). For n = 1.5, more than 
80% emission energy can be collected within 9 ~ 11°. 
For n = 1.9, very good directionality can be obtained, 
with over 50% collection efficiency for < 0.7°. 

For widely used microcavities with the refractive index 
not equal to 1.9 (e.g., Silica, Calcium Fluoride, Lithium 
Niobate), the optimal emission directionality can be ob- 
tained by optimizing the shape of the microcavity (e.g., 
using deformed cavities) and the position of the scattcrcr 
(e.g. embedding the scatterer inside the microcavity). 
Note that this optimal refractive index is irrelevant with 
the size, shape and material of the scatterer. 



FIG. 7: (Color online) (a) Half energy angle Q\/2 vs. n for 
4> = 0°, 90°. 
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V. CONCLUSION 

In summary, based on cavity QED approach, we an- 
alytically investigate the microsphere-scatterer coupling 
system in which the high-Q WGMs can be efficiently ex- 
cited through free space and the resulting laser is capa- 
ble of emitting with high directionality. In this system, a 
subwavelength scatterer placed in the vicinity of the mi- 
crosphere serves as an interface between the input light, 
the WGMs and the output light. We take advantage of 
the microsphere itself to focus the input beam with small 
spot area, and collimate the output beam with ultra- 
small divergence angle. Our results show that the high- 
Q WGMs can be excited with efficiency larger than 10%. 
More importantly, the half-energy angle of the output 
light can be as narrow as 0.7°, which is a great improve- 
ment over the 2D microcavity lasers [HI, [33] ■ This holds 
great potential for novel micro-sized laser sources and has 
broad applications in Micro/Nano photonics. 



FIG. 8: (Color online) Energy density p(0) and energy ratio 
P(0) as functions of different n and <j>. (a)-(b): n — 1.5; 
(d)-(e): n = 1.9. 

In Fig. [7] we plot the half-energy angle as a function of 
the refractive index for <fi = 0° and 90°. Note that these 
two cases set the lower and upper bounds, as inferred 
from the emission pattern in Fig. EJa). Remarkably, 
0i/2 can be less than 1° (for refractive index around 1.9), 
which indicates much better directionality than previous 
predictions based on other mechanisms, to the best of 
our knowledge. 
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Appendix A: Derivation of Quantum Langevin 
equations 

Starting from the total Hamiltonian H\ ot (Eq. JT])), 
the Heisenberg equations of motion can be derived as 



8 



+ 00 

= -iuia m - J g in (u)b(u})duj + ig x a m > - i^gijCj, (Al) 



dt 

db(u) 
dt 

d 
dt 



j 



= -iub(u) + g- ln (w)y j a m , (A2) 

j n 



I 

The eigenmodes of the system are the symmetric and accw)/ Then the above equations can be rewritten 
antisymmetric standing modes, given by a± = (acw ± as 



+00 

da f 

— = —iu)ia + — \/2 I g- m (uj)b(u))dLU + i2g\a + — i\f2 9i.j c j h (A4) 



da- 
~dt 



db(oj) r- 
— — — = —vxibyuj) + V2gi n {(jj)a 



dc^ 
dt 



Formal integrations of Eq. (|A6[) and (|A7j) yield 



\ 

_l 



(A5) 
(A6) 

-iV2g hj a + . (A7) 



b{uj) = e- lul{t - to) b (uj) + V2. 9in (cj) j e-^-^a+^dt', (A8) 

t 

Cj = e-^ ( *- to) Cj -o - iV2g hj [ e-^^-^a+^dt', (A9) 



I 



where bo(u), Cj,o denotes the value of b(uj), Cj at t = to, arise from the interaction with the WGMs. 
respectively. In both equations the first terms represent Substituting the solutions into Eq. (| A4|l . wc finally 

the free evolution of the modes while the second terms obtain 



= - 2g l )a + - (K in + kr)o+ - V2K~b- m - £, (A10) 



I 

where represents the input- WGMs energy coupling strength, 

*- = - 25l) (A11) V = J= r bo^e-^d. (A12) 

V27T J-00 
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describes the input field, 

KR = 2-^//-.^ - cjj + 2 5l ) (A13) 

denotes the damping of the WGMs induced by the scat- 
tering to the reservoir, 

^v^V.^e-^'-'^.o (A14) 

is the noise operator relate to the reservoir. In deriving 
Eq. (|A10[) . we have used the Markov approximation [46| . 



Then the transmission can be obtained as T = 
(bl ut b ou t) / (bj n b i n ) . For optical frequency, room tempera- 
ture, the initial states of the reservoir modes are almost 
all vacuum states. Thus the expectation values of the 
noise operators can be neglected. Therefore, we obtain 
the minimum value of the transmission 

T ^=( lT7:~lT ) 2 - (B3) 



Appendix B: Derivation of excitation efficiency Thug the excitation efficiency r, = 1 - T min can be ob- 

tained as 

Taking the intrinsic decay rate kq of the WGMs into 
account, and using the input-output relation (4(| H3] 

&out = kin + V2~K~a+, (Bl) 4K in (^ + 2k r ) 

(k + 2« in + 2k r ) 2 ' 1 ' 

we obtain 

, , 2K; n &i n + \/2Ki n t; ,„„, 

0out=0 in + T7 — : 77T— ; r. (B2) 

i(u - CJi + 2g 1 ) - {f + K in + kr) 
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